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Abstract Consistent with the predictions of Kibble and Zurek, scaling behaviour 
has been seen in the production of fluxoids during temperature quenches of super- 
conducting rings. However, deviations from the canonical behaviour arise because 
of finite-size effects and stray external fields. 

Technical developments, including laser heating and the use of long Josephson 
tunnel junctions, have improved the quality of data that can be obtained. With new 
experiments in mind we perform large-scale 3D simulations of quenches of small, 
thin rings of various geometries with fully dynamical electromagnetic fields, at 
nonzero externally applied magnetic flux. We find that the outcomes are, in prac- 
tice, indistinguishable from those of much simpler Gaussian analytical approxima- 
tions in which the rings are treated as one-dimensional systems and the magnetic 
field fluctuation-free. 

PACS numbers: 74.40.Gh, 05.10.Gg, 98.80.Bp 



1 Introduction 

Many systems of interest in condensed matter physics exhibit a second-order 
phase transition that can be studied in the laboratory. For such transitions one 
would expect, under adiabatic change, that the correlation length of the system 
could grow to be as large as the system size. The system therefore orders itself si- 
multaneously and in the same manner everywhere. In practice, however, causality 
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(the finite speed at which order can propagate) restricts the size of the domains 
within which the system can order itself and the correlation lengths are rather 
smaller. The faster the transition, the shorter the largest correlation length that is 
reached, and hence the smaller the domains. 

If the phase transition in question results in a spontaneously broken symmetry, 
different choices of the order parameter can be made by domains that have been 
causally disconnected while the system was out of equilibrium. If there is frus- 
tration then topological defects will form, which may persist long enough to be 
detected. 

Such spontaneous formation of topological defects from phase transitions out 
of equilibrium is frequently termed the Kibble-Zurek (KZ) scenario. The process 
outlined above was discussed by Kibble-^ in the context of Grand Unified The- 
ories (GUTs) describing the physics of the early universe. Zurek independently 
suggested that the same process could occur in condensed matter systems^. 

The most straightforward (and elegant) manifestation of the KZ scenario in- 
volves a linear temperature quench from one phase to another. Let us parameterise 
this linear quench by the time taken Tq = T c (dT / dt)j l , where T c is the critical 
temperature of the second-order phase transition in question. 

If the system is initially homogeneous and isotropic it is proposed that the 
initial defect separation £, will be given by a scaling form 



where and To are system-dependent. The exponent a, however, can be de- 
rived from critical exponents valid for the adiabatic transition and so the essen- 
tial university of the phase transition can be seen even in this out-of-equilibrium 
behaviour. 

Unfortunately, we have yet to see evidence for defects in cosmological obser- 
vations, although they are predicted in all reasonable models that incorporate uni- 
fication. However, Eq. (fl3 has been tested in a wide variety of different condensed 
matter systems. Restricting our attention here to experiments involving type-II su- 
perconductors, it has been demonstrated 5 '^ that such arguments are applicable to 
Josephson tunnel junctions and to isolated superconducting rings 7 ; these experi- 
mental results are supported by simulations^. 

There is one significant difference between cosmological and condensed mat- 
ter systems; whereas cosmological phase transitions take place in effectively in- 
finite volume, those in condensed matter experiments are limited to finite - and 
often quite small - sizes. This can lead to edge effects, and allowances should be 
made for the different physics close to the boundaries of such systems. 

We shall concentrate here on fluxoid formation experiments in rings of Nio- 
bium. If the ring is sufficiently narrow, then it is plausible that the ring behaves 
like a one-dimensional object with finite size and periodicity. This accords with 
the KZ picture, for which Eq. (Q} is now understood as a perimeter law. A pre- 
vious study tested these assumptions and examined how the KZ scenario unfolds 
for slow quenches in small ID systems with periodicity—. 

Here, this work is extended in two principal ways. 

First, we extend our simulations from one dimension to three dimensions, to 
describing a flat annulus of superconductor in three-dimensional space, in which 




(1) 
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the electromagnetic field is fully dynamical (impossible in one or two-dimensional 
simulations). This allows us to test different ring geometries (the KZ picture im- 
plies only perimeter scaling behaviour, indifferent to ring shape and width), and 
the breakdown of the KZ scaling laws for small systems. In addition, with dynamic 
electromagnetic fields, we can determine the likelihood of observing Hindmarsh- 
Rajantie magnetic field freeze-out 1 1 , which also goes beyond the KZ scenario. 

Secondly, we discuss the effect of an external magnetic field on the supercon- 
ducting phase transition and compare with simple theoretical estimates for how 
this biases the fluxoid trapping rate. Experimentally, it is impossible to eliminate 
stray magnetic fields completely; if we can understand the dependence of fluxoid 
production on the applied field than we can reduce errors due to stray field bias. 



2 Analytic approximations 

As we say, realistic systems can give deviations from the simple behaviour of 
Eq. (QJ because of finite size and stray fields. To provide an estimate of the likeli- 
hood of finding spontaneously produced fluxoids within our ring taking these into 
account, against which we can test our simulations, we adopt simple Gaussian 
ansatze, assuming that we can treat the annulus as a ID system. 

Let us denote the winding number of the complex order parameter field <P 
(which mimics Cooper pairs) along the ring by n. That is, if the ring has circum- 
ference C, the total winding number is 



where 8 is the change in phase around the ring. The probability of seeing n flux- 
oids or anti-fluxoids is /„ = p n + /?_„. We shall concentrate on the slow quench 
regime in which f\ <C fi and can ignore the production of more than one unit of 
flux. 



2. 1 Finite size effects 

We have discussed this elsewhere 10 but, in brief, we assume that the winding 
number density is a Gaussian random variable. If C is the circumference of the ring 
the relevant quantity is S, /C. This determines whether the ring is 'large' or 'small' 
according as the ratio is small or large. Finite size effects arise either because the 
ring is physically small or because the quench is slow. 

In the absence of an external field this Gaussian ansatz is sufficient to repro- 
duce the KZ scenario for large rings, with exponent a = 0.25. However, for small 
rings we find exponential damping, in which Eq. (Q} is replaced by 




(2) 



ln/i 



A| 2 /C 2 + const. 



(3) 



The coefficient A is not simply estimated. Eq. (O is in contrast to the proposal that 
the slope be doubled 7 , which may arise for a disc with a small central hole, but 
which we have not tested here. 
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2.2 External field effects 

Now make the slightly different assumption that Arg <P itself has Gaussian correla- 
tions around the ring. In the presence of a field, it is proposed that the accumulated 
phase change 8 then has a Gaussian probability distribution G(0;0f) 

» -p-2# (4) 

V27TCT 2 2(7 Z 

where we have assumed that the mean is biased from zero by the applied nor- 
malised flux </>f through the ring; the deviation a depends on geometry and quench 
time. With such a distribution, the probability p„ of ending up with a given wind- 
ing number n is given by 

/Tl+lnTZ 
dd G(0;0f). (5) 
-7I+2MZ 

In the presence of a stray residual field, </> r , we replace </>f by 0f — r . The probability 
of trapping n fiuxoids can be found to be 



Pn(<k) = 2 



„0f- </>,-- n + 0.5 . </> f -(/»!■- n - 0.5 
erf erf 



(6) 



where the dependence on quench time and geometry is parameterised by s. Sim- 
ilar expressions can be derived for the probability of trapping other numbers of 
fiuxoids. This calculation unfortunately does not allow us to determine the depen- 
dence of s on Tq. Nevertheless, we can fit data for different applied magnetic fields 
at fixed Tq to Eq. $6$ with only s and </> r as free parameters to yield /q. We will test 
the validity of this ansatz (with </> r = 0) in our simulations. 



3 Simulations 

To test the approximations outlined in the previous section, and investigate the 
extent to which the physics of a superconducting ring is one-dimensional, we car- 
ried out simulations. In our system, illustrated in Fig.Q] the superconducting ring 
(in which the complex scalar field <P is defined) is surrounded on all sides by a 
box in which the standard U(l) vacuum electrodynamics is simulated. Although 
the superconductor is thereby taken to be a thin planar film, the three-dimensional 
simulation allows for any nontrivial correlations of the magnetic field to play a 
role in addition to the KZ mechanism. We can therefore study both mechanisms 
that may be thought to affect fiuxoid formation in this system. 

We use a gauge-invariant Langevin formulation for gauge and scalar fields, 
making use of Krasnitz's work on formulating Langevin equations for systems 
with first-class constraints^. 

A noncompact formulation of the U(l) gauge field is used, so that it can take 
any value on each link rather than being defined up to 2n. The Hamiltonian for 



Fig. 1 Simulation setup. In lattice units (we typically take a = 0.5), the ring - the shaded region 
marked &?, - is always of a constant width of 5. The border around the ring is also 5 units thick; 
the size Ax x Ay of the central void therefore determines L x and L v ; we take L z = 5. The external 
magnetic field 5 e xt is applied in the direction shown. Periodic boundary conditions are used in 
all three directions. 



the system is 

X,l 



jk 



•2 £ Re<P*(x)Ui(x)<P(x + i) 



n*n 



<p*<p + ^(<p*<p) 2 



(7) 



where t/,(x) = e iaeA >( x \ and the finite differences are 



Affix) 



f(x + ai)-f{x)_ 



Arf{x) = f ^- f{x - at) . (8) 



Note that <P is only defined on a 2-dimensional superconducting ring region de- 
noted M. The only gauge fixing that has been carried out is Ao = 0; 'temporal 
gauge' . This yields the Gauss constraint 



2e, 



£zT£, = — Im <P*n. 



(9) 



The generators for the noise have been chosen to be the gauge-invariant observable 
quantities {E,, l^l 2 }, in line with previous work on superconducting films that 
made use of the Krasnitz metho d 13 ' 14 . This physically-motivated choice is not 
unique. 
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Defining covariant derivatives 



the stochastic field equations are then 



(10) 



Ai^Ei+tfEi+r) (11) 

£ = £ e ijk e klm AjA+A m {x) - -Y Tm&*(x)Df&(x) (12) 

4> = n (13) 

J7 = £Z)rD+*(jr)-m 2 <f'*-A(***)*-*(j3 /T ^|*| 2 +r n ) (14) 



where F and F n are the Gaussian noise terms for the gauge and scalar fields, 
respectively, satisfying fluctuation-dissipation relations 

(r(x,t)r(x!,t')) =2l3T8 {3) (x-x')8(t-t') (15) 
(r n (x,t)r n (x'/)) =2p n T5 {2) (x-x')8(t-t'). (16) 

The dimensionalities of the delta functions differ due to the scalar field lying in a 
plane. We will take j3 = 0.5, f5 n = 0.25, T = 0.01, m 2 = 1, e = 0.1 and k = 0.1; 
this places us in the type II regime with K = w 6.3. We have checked that 

our results do not depend on lattice spacing, but for convenience we generally 
work with a = 0.5. The values of j3 and (5 n were chosen such that the system 
behaves in a heavily damped manner after the quench. 



3.1 Boundary conditions and External field 



Periodic boundary conditions are used at all times. This forces the total magnetic 
flux through each of the boundaries to be zero. Where we wish to impose an ex- 
ternal field, we do so by twisting a single plaquette located at an arbitrary position 
xq in the x — y (1 — 2) plane of the gauge part of the Hamiltonian— . The effect of 
this is to add the term 



H ext = aY,(A l (x)+A 2 (x + al)-Ai(x + a2)-A 2 (x))(l) ext S i2) (x-XQ) (17) 

X 



to the Hamiltonian, where </> ex t is the applied flux. We are free to make this greater 
than 2n, because we have used the noncompact formulation of the gauge field. 
The equations of motion must of course be modified appropriately. It is possible 
to think of this external field as a Dirac string running through one plaquette on 
each slice in the z-direction. 
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3.2 Quench protocol 

The stochastic equations given above are used to simulate the system throughout 
(although the same results could be obtained at less computational cost by using 
an initial thermal ensemble obtained from a Monte Carlo simulation and then 
evolving with overdamped dynamics). Given an initial 'cold' system, the system 
was thermalised until the magnetic flux density in the z direction was constant 
throughout the lattice and the time-averaged total energy was no longer varying. 
A linear quench was then carried out to take the system into the superconducting 
phase. Starting from t = — Tq, m 2 was changed from its initial value as follows 



At t — 8tq the simulation was stopped and the winding number n was measured; 
the total magnetic flux trapped in the ring was also measured (as a cross check). 



The system was simulated as described in Section [3j Each data point is the result 
of 600-1000 individually thermalised runs; the results presented in this paper rep- 
resent about five thousand hours of computer time. Error estimates were obtained 
from a bootstrap resampling of the results. We verified that the results did not 
change when the border around the film was enlarged, nor when L z was increased. 



4. 1 Finite size 

With S ext = 0, we first sought to see if the exponential falloff in trapping probabil- 
ity predicted 10 by ID simulations in Eq. OJ persisted for the more realistic theory. 
It does. In Fig. [2] we show the crossover from KZ behaviour, indistinguishable 
from the falloff obtained previously 10 . Differing aspect ratios were tried for the 
central hole for the same perimeter, ranging from square to an elongated 4: 1 rect- 
angle, but no deviation from the perimeter law behaviour of Eq. (QJ was detected. 
In this parameter range, then, we expect that the KZ mechanism (and its extension 
to small volumes) will be the only detectable influence on fluxoid formation, with 
no discernible evidence for magnetic field correlations. 



4.2 External fields 

The results for nonzero external field are shown in Fig. [3] with fits of the form 
in Eq. (|6]l also shown. Despite the simple assumptions going into this ansatz, it 
agrees with the results of the simulations very well indeed. We tested the ansatz at 
various Tq and found that it worked well throughout the range shown in Fig. [2] 




-Tq < t < Tq 
t > Tq 



(18) 



4 Results 
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Fig. 2 Probability of observing one defect, /|, as a function of quench time Tq. Results are 
shown for two aspect ratios. As observed in the one-dimensional simulations, a change in be- 
haviour from KZ scaling to exponential suppression occurs. The power law fit shown is to the 
first four points only and yields an exponent a = 0.26 ±0.05. Incorporating additional points 
produced a poorer fit; we associate these with the exponential damping regime discussed in the 
text. 



5 Discussion 

We have simulated the effect of a temperature quench on a two-dimensional su- 
perconductor with its central area removed, embedded in three-dimensional space 
(see Fig. QJ as it is driven through its phase transition. The results of our simula- 
tions are two-fold. 

Firstly, in the absence of an external magnetic bias field, for the trapping rates 
studied there is spontaneous production of fluxoids which, for rapid quenches, 
accords with the KZ scaling behaviour of Eq. (QJ and which, for slow quenches, 
shows exponential falloff; there is no perceptible influence of the magnetic field 
freezing out. Recent experiments have also struggled to show this effect 16 . Fur- 
ther, the fluxoid density depends only on the inner perimeter and not on the as- 
pect ratio of the superconducting rectangle. In all of these regards, there is no 
difference between the full theory and what we have observed in an idealised one- 
dimensional ring 10 with Gaussian winding number density. 

Secondly, experimentally it is difficult to eliminate stray magnetic fields when 
carrying out temperature quenches of superconducting rings. Such fields artifi- 
cially increase the likelihood of seeing fluxoids. Although the simulation of Fig. [3] 
is for zero stray residual field </> r we now know how to proceed to eliminate their 
possible effects. We apply external fields and identify the minimum (true) value 
of f\ by matching the profiles to those of Fig.[3j Again, the full simulation is well 



9 




Fig. 3 Probability of observing zero (po), one (pi) and two (P2) defects as a function of nor- 
malised external magnetic flux (j>f, with Tq = 120. The curve is a one parameter fit of the form 
given in Eq. ® for p\, and yields s = 0.39 ±0.01. 



represented by a theoretical Gaussian ansatz that models the effect of an external 
(non-dynamical) field on fluxoid production. 



The relevance of this analysis is that experiments are being prepared to provide 
a rigorous test of the simple KZ scenario for superconducting loops, originally 
proposed by Zurek more than 25 years ago but, as yet, only partially attempted—. 
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